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In this paper, we investigate the perturbations in matter bounce induced from Lee- Wick lagrangian 
with the involvement of non-minimal coupling to the Einstein Gravity. We find that this extra non- 
minimal coupling term can cause a red-tilt on the primordial metric perturbation at extremely large 
scales. It can also lead to large enhancement of reheating of the normal field particles compared to 
the usual minimal coupling models. 

O ! PACS numbers: 

^ \ I. INTRODUCTION 

>; 

^ . Non-singular bouncing cosmology [15 has become one of the important alternative scenarios of inflation [2[ as a 

description of the early universe. In terms of a bounce scenario, the universe travels from a contracting phase to an 

^v ,, expanding phase through a non-vanishing minimal size, avoiding the singularity problem which plagues the Standard 

psj ' Big Bang theory [3] or inflationary theory |^]. Moreover, a bounce can occur much below the Planck scale, eliminating 

\ the "Transplanckian" problem of which the wavelength of the fluctuation mode we see today will be even smaller than 

'^^' • Planck scale and thus in "the zone of ignorance" where high energy effects are robust and Einstein equations might 

^w ^ be invalid |5|. In bounce models, fluctuations were generated in contracting phase and transferred into expanding 

\^ . phase through bouncing point, which can give rise to a scale invariant power spectrum as expected by the current 

observational data 6]. In order to connect the perturbations in contracting phase and expanding phase at the crossing 

O point, the joint condition of Hwang- Vishniac 7] (Deruelle-Mukhanov Q) can be applied. Besides, bouncing scenario 

I ■ can also be originated from non-relativistic gravity theories [9| and give rise to gravitational waves with running 

H \ signature and large non-Gaussianities |ig |. 

*T* ' In the previous studies, a kind of "matter bounce" has been proposed where before and after the bounce, the universe 

j^ \ could behave as non-relativistic matter. This scenario can be realized by a kind of "Lee- Wick" (LW) Lagrangian which 

~~^ • contains higher derivatives of the field, or equivalently by two scalar fields with an opposite sign of potential. This 



^ 



^ 



> 



>< 



S 



^ ' scenario can give rise to scale invariant power spectrum at large scales, which is consistent with the data [ll|. However 



one may expect that there are some non-minimal couplings of the matter in the universe to Einstein Gravity in the 
early universe [33]. In fact, it is argued that the existence of the non-minimal coupling term is required by the quantum 
[~^ , corrections and renormalization 12] in quantum gravity in curved space-time. This effect may, by means of modifying 
lO ' Einstein equations, alter the previous result and leave some signals on the power spectrum that can be detected by 
Cn _ observations such as CMB analysis. The non-minimal coupling can also be reduced from higher dimensional theories 
t^^ . such as brane theory and can get rid of the Big Bang Singularity, which lead to a bouncing universe [1^ i or lead to 
f^ ' the cosmic acceleration, which can be utilized as infiaton in the early time [19| and dark energy at current epoch |20| . 
^D , This paper aims at investigating the perturbations in the matter bounce model involving the non-minimal coupling. 

^~~^ ' The paper is organized as follows: we first review the model of matter bounce in Sec. II. After that, in Sec. HI we 
^ , take the non- minimal coupling term into account. We investigate the perturbation through the process in detail, and 
show the solutions for each stage. We also analyze the particle production due to the resonance. All the analytical 
calculations are done with help of numerical computations. Finally Sec. IV contains conclusions and discussions. 



II. REVIEW OF THE MATTER BOUNCE MODEL 



In this section, we will start with the matter bounce model carried out in [11| . This model consists of a single scalar 
field with a high derivative term of a minus sign. Due to this term, the equation of state (EoS) w is possible to go down 
below —1 and violates the null energy condition, which behaves as a "Quintom" matter [2lj and makes it possible to 
realize a bounce. It is useful to have such a term added in the lagrangian. For example, in the well known Lee- Wick 
theory which was constructed by T.D. Lee and G. C. Wick ^^ (see also [23] for the extensive "Lee- Wick standard 
model" ) , the higher derivative term is used to cancel the quadratic divergence of the Higgs mass and thus address the 
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"hierarchy problem" . It is also possible to construct an ultraviolet complete theory which preserves Lorentz invariance 



and unitarity in terms of Lee- Wick theory [24 1 . 



We begin with the Lagrangian in our model to be in the form: 

£ = ld,$d^$ - ^id'4>r ~ l^^^ - Vi$) , (1) 

where m is the mass of the scalar field (f) and V{(j)) is the potential. A higher derivative term with minus sign is 
introduced with some new mass scale M. For the Higgs sector of Lee- Wick theory, the hierarchy problem is solved if 
we require that M ^ m. After some field rotations, one can write down the effective Lagrangian: 



^d,,<j>d'^cj, - ^d^4>d''4> + \m^4>^ - \^ 



C = -d^^cjid^cj) - -^^^^^'4> + -M^^^ - -m^cf)'^ - ¥{(/), 4>) , (2) 



where (j> is some auxiliary field and (f> is defined as = (p + (j) |34l |. Here the mass matrix of the fields has been 
diagonalized due to the rotation. Usually there may be some interaction terms between the two fields, or some higher 
order self-interaction terms like (/)^, 0*, cfPcj)'^ and so on, but here for simplicity and without losing generality, we will 
have all of them ignored by setting V{(j), (p) = 0. In framework of Friedmann-Robertson- Walker (FRW) metric. 



ds^ = df - a{tfdy? , (3) 



it is easy to get the Friedmann equation as: 



H' = ^[\^' -\4^ + \m^4>' -\m^'] , (4) 

where k^ = SttG, and the equations of motion of the two fields are 

4> + 3H,j) + m'^(l) = , 

-I- 3H^ + M^4> = , (5) 

respectively. 

Let us now take a close look at how the model works to give rise to a bouncing scenario. The bounce happens by 
the conditions H — and H > 0, which requires that the total energy density vanish at some point in the universe 
evolution. Starting off in the contracting phase, both the two fields and 4> oscillate around the extrema of their 
potential with a growing amplitude due to the anti-friction of the negative Hubble parameter. Therefore, the energy 
densities of both fields will grow as a^^{t), which behaves like non-relativistic matter. Note that the energy density 
of (j> is negative and its growth will cancel that of the total energy density of the universe. At the beginning, we 
assume that it is subdominant than that of </), but as M ^ m, it will grow faster and eventually catch up, canceling 
all the positive energy density with its kinetic energy overwhelming that of cj), which makes H > and cause bounce 
to happen. Thus we can see that before and after the bounce, the universe can be viewed as "matter domination" 
with the effective equation of state w — 0, while in the neighborhood of the bounce point, the equation of state goes 
down to — oo |35| . It is a good approximation to parameterize w in such a way for the following calculations. 

In the following section, we will reinvestigate the scenario by taking into account the coupling of the field to Einstein 
Gravity. We assume that the normal scalar couples to the curvature through some non-minimal coupling term such 
as i^i?(/)^. We will analyze in detail the evolution of the perturbations in presence of this term, including the metric 
perturbations as well as the effects on particle production. To support our analytical calculations, we will also perform 
the numerical investigations. 

III. ADDING NON-MINIMAL COUPLING TERM TO THE LW BOUNCE 

A. background 

In this section, we investigate the system with the inclusion of the non-minimally coupled term. Generally, we can 
modify Lagrangian (H)) in a very general form such as: 



l^^^^'^-^^i^'^r-k 



C = -d^^d^^ - TTTT^id^f - -m202 _ v{^) - f{R)h{(f>) , (6) 
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where f{R) and h{(j)) are functions of the Ricci scalar R and scalar field (j), respectively. By writing it in terms of i 
and (j) in analogy to ^, we can always have: 



Id.^d^^c^ ~ ld,4>d^^ + ^M^^ - i. 



C = -d,^d^^4> ~ 0^^d^^+-M'^^ - -m'ct,^ - y(0,0) - h{R)hm - /2(i?)/i2(0) , (7) 









+ 3H0 + m^<\> + t,R(i> 
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where fi{R), f2{R), hi{(f>) and h2{(j>) are inherited from f{R)h{(f>) in the original Lagrangian. Usually, these functions 
are very complicated and difficult to analyze. For the sake of simplicity, we adopt a kind of parametrization form in 
which the coupling of (j> has a very simple form such as ^S^RcjP with a coefficient ^, while that of 4> vanishes. The 
parameterized Lagrangian is: 

From this Lagrangian, we obtain the Einstein Equations as: 

j^G^i- =5mi'[o'?^^-o<^ +-m202--M2(/,2]_v^<?;.V^<?;.+V^0V^0-2^5^^(,^2_,_^Q^)_^2^(V^0V^0+0V^V^0) , 

(9) 
and the equations of motion of </) and </>: 

+ 3H0 + m^^ + t,R<^ = , 

(10) 

2ij2). 

(11) 



(12) 

Generally, the additional term will affect the background evolution by modifying the field equations. This may 
make the analysis sophisticated and undoable. Moreover, if the modification is too large, the bounce may even be 
prevented. However, for the case of weak coupling, it is still suitable for us to preserve the bounce and take assumption 
that the background evolution is hardly affected. Thus it will be safe to use all the results obtained in the above 
section without losing fidelity. 

To be more accurate, we test the background evolution numerically. We choose the parameters and initial conditions 
to be the same as in [ll|. We find that up to |^| ~ 0.05, the background evolution will not be affected very much 
and a bounce can still happen. But there may be a little difference, where the universe may enter into a small stage 
with w ~ — 1 before bounce, as have been shown in Fig. [TJ This is easy to understand from the field equation. In 
the additional term, ^_R can be viewed as an effective correction to the original mass squared vr? . At the initial stage 
where R is not very large, the correction is negligible, and the field evolves as it does without the additional term. 
But as the universe evolves like pressureless matter, R will get larger and larger. When the Hubble parameter is 
about to its maximal value, H ^ m^ the second scalar will begin to dominate the universe where H transfers from 
negative to positive. At this moment, from the definition of i?, we can see that it suddenly jumps by the amount 
of 2|ij|. This behavior makes the energy density dominate the universe again for a short while, leading to ii; ~ — 1 
in this period. As the energy density of cj) grows up and totally dominates the universe, the equation of state goes 
down to be below —1 and drives the bounce to happen. At the bouncing point where the total energy density of the 
universe vanishes, one can see from Fig. [1] that the equation of state goes down to negative infinity, and as stated 
before (in the footnote), it is not a real divergence. After the bounce, the EoS will come up to above —1 again as the 
time reversal of the pre-bounce process. 

B. perturbation 

In this section we will focus on the perturbation evolution of our model. We begin with the perturbed metric, 

ds^ = a(r/)2[(l + 2<i>)dr]^ - (1 - 2*)dx2] , (13) 
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FIG. 1: (Color online) Evolution of the background equation of state parameter w and scale factor a{t) in a non-minimal coupling 
matter bounce model as a function of cosmic time t (horizontal axis) . The background fields are plotted in dimensionless units 
by normalizing by the mass Mrec = W~^mpi, and the time axis is displayed in units of M~^].. The mass parameters m and 



M are chosen to be m = b.OMrec and M = lO.OMrec- The initial conditions are 



-370.5UMre 



-6752.791M4 



-4.818 X 10""Mr 



-8.837 X 10~"M;, 



in the conformal form. Here rj = f — is the conformal time, and <J> and ^ are scalar perturbations of the metric. Note 
that with existence of the non-minimal coupling term, $ and ^ are no longer equivalent. Similarly, the fields can be 
represented as homogeneous and fluctuation parts as: 



(14) 



In the following, we will neglect the subscript "0" and take and (j) as background components. Then we will solve 
the perturbed Einstein Equations to see their evolutions. 



1. general solution 



From the perturbed Einstein Equation, 
we obtain the following equations for perturbation variables 



SG^,^ = -K ST^i, , 



$ = * 



F K-2 

^ + (H )$ = _(- 



6F 
~F ' 



5F HSF 
2F^ 2F 



(15) 

(16) 
(17) 



^ + 3(i/ + ^)^ + I y [(0' + 12^iJ#) + 2^^'] - {2H + m^) - 2^ 



5F 3 



{2H + iH'^ 



+ ( 



F 12 






F 



p )5<P~-q + 2HS,cj,)5<l,, 



K^m^d? F 

^ * 

F \ 




6F F6F 



2F 2F2 



(18) 



,k 



Scf) + SHScf) + \ (-) -f m^ -f- 6^{H + 2H^) 



6$<^* + 24^iJ0* + 30* + Q^H(j)<i> + 0$ 



4Cfc20* 



2^(P[6{H + 2H^) - {-f]<^ - 2{m^(t) + ^R(j))^ , (19) 
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where we have defined F = 1 — K^^cfP' . 

Since the equations are too sophisticated and there is no hope to solve it directly, we resort to the Einstein frame 
by a conformal transformation of the metric [231 ■ 

9^J,u = Fg^, . (20) 

The action ([8]) is then transformed into the form: 



-2^2 2F^ ""^ ^ 2F ''^ "^ F2'2 "^ 2 

We introduce two new variables, x and x, which are defined as: 

Vi - (1 - eo^^^.^^ 
X ^ ■ 



X = I id4>, (22) 



F 

3.6 . 
F 

and rewrite the above action as: 

5 = 1 d'^V~9[^ + \d,xd''x - l^xd'-X + ^(iM^ixf - Im'Hxf)] ■ (23) 

The perturbed metric will be: 

ds^ = a'^{r])[{l + 2^)drf - (1 - 2#)dx2] , (24) 



where a(rj) — \'Fa{r]). 

For the perturbed metric in Einstein Frame where the off-diagonal components of the perturbed Einstein equations 
vanishes, we have $ = ^. This can also be obtained when we consider up to the first order approximation of the 
perturbed metric, which gives $ = <& + |^, and ^ = \l/ — |^ where relation between $ and ^P has already been given 
in Eq. p^ . Moreover, Eq. p7)) gives the direct relationship between $ and the field perturbation, 

... !,.!_ 2FHaVF^y ^^^^ 



^ ^ ^ ^ ^2„[1 _ (1 _ 6^)K2e02] ' 

which is convenient for us to calculate $ in Einstein frame and then transform into the original Jordan frame. 
The equation for <t is: 

I*" + 2{n - ^)$' + 2{n' - n^)^ + fc^d = K^{2n + ^)x'Sx , (26) 

A A A 

where H = a' /a. The right-hand side (r.h.s.) of the equation is assumed to be small and negligible except at the 
bounce point. Moreover, it is also convenient to define the curvature perturbation ( in terms of <&: 

c^^+ ^J^ ^ (^'+m, (27) 

which is expected to be conserved on super- Hubble scales in the inflationary universe by the condition that the entropy 
perturbations are small and (1 + w) 7^ 0: 

{l + w)C'=0. (28) 

At the next step we will solve Eq. (l26l) for each stage to get the perturbations $, 'J and Scf), and sec how it is 
effected by the non-minimal coupling term. 
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2. contracting phase 

In contracting phase, the universe is dominated by the normal scalar (p. Since F deos not vary significantly, we 
roughly have x ~ which does not play an important role in the evolution. Then by neglecting the r.h.s. of Eq. (j26l) . 
we consider it as a homogeneous equation. It is useful to define a new variable u = a^/x' ^^^d rewrite Eq. (|26|) as: 

u" + ik^~^)u^O , (29) 

a 

where 9 = 'H/{dx')- The solution can be split into two limits: short-wavelength (fc^ ^ 9" /9) and long-wavelength 
(/c^ ^ 9" /9) where 9" /9 is the potential. For the short-wavelength perturbations where the potential term can be 
neglected, we have: 

u = Ce^'^" + c.c. , (30) 

where C is a constant and c.c. stands for the complex conjugate. For long-wavelength perturbations, one gets: 



u = 9{S + DJ'^). (31) 



From the above argument of $ and 4* as well as the definition of u, one can find the solutions for perturbation variables 
in terms of u which are: 

*^is-...j;^4)'. (32) 



a 



2aF{n^ -W) a9 



uH F' uH^, 

S^^ 2aF{n^-n') ~&9 



^ - ^+ ..^...o .,J —y^ (33) 



C^y , (34) 



aCW^-n') a9 

where in the last formula we have neglected the contribution arising from 64> in (|25p . and we have used the approx- 
imation x'^ - 2(7^2 - 'H.')/k^. Substituting ^ or ^^ to ((221), ^ and dM]), we get these variables in Jordan 
frame: 

'J>-C^\/5^ -^-i^(^+-^ + f)K--> (35) 

H [ 2H 

^-C^I-.^\-^+liH+A + j)}+c.c., (36) 

-H 2H 



OO 



S^ = ^^Vi^(^+^ + f)+-- (37) 

y-H 2H 

for the short-wavelength case, where H is the Hubble parameter in Einstein Frame and the symbol "o" denotes cosmic 
time derivative in Einstein frame (for details see Appendix A), or 

6(j) = -^{S^^D f aFdt) (40) 

for the long-wavelength case. The sketch plot of fiuctuation modes compared to 9" /9 are shown in Fig. [H 
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FIG. 2: (Color online) A sketch of the evolution of fluctuation modes compared to the potential 6" /9 in a bouncing universe. 
The vertical axis is co-moving spatial coordinate, and the horizontal axis is the conformal time. The orange line denotes the 
fluctuations with small comoving wavenumber k while the blue line is for those with large k. 

For the case that the universe evolves with some constant equation of state w, the scale factor can be parameterized 



a{t) - tP 



(41) 



where p is defined as p = 2/[3(l + w)]. This will be the case in the whole process except for the bouncing point. 
During the contracting phase before the non-minimal term dominates the universe, the fields oscillate around their 
extrema which made the universe behave like non-relativistic matter. The average vakie of the equation of state of 
the universe is < w >= (see from Fig. [ij, and the scale factor a{t) scales as t^/^. Since in this period the factor F 
does not evolve significantly, it can be approximately viewed as a constant. Therefore the scaling of $, ^ and (50 are 
given as: 



$ 



2C 



4 1 _.ia 

e "0 

1 ik 



-)e- 



-fs 



c.c. 



for short-wavelength case, and 



$ 



* 



25'; 



6D- 



3aoFt^ 5k^ 



o^ 



qqF V 3k 



:r^ + 



6D, 



4 
3^ 



(42) 
(43) 

(44) 
(45) 



for long-wavelength case, where oq is some normalization constant. Here the subscript "c" denotes the contracting 
phase and superscript " — " means the matter-like region which happens earlier. Note that $, ^ and S(j) have 
magnitudes of the same order which should be much less than rupi , so the difference between $ and ^ (proportional 
to 6(j)) will be severely suppressed. 

However, when the non- minimal term begins to dominate the universe, the equation of state approaches to —1. 
During this period one could define a slow roll parameter e = —H/H^ which should be very small. Thus we 

can have a — ait^''^ and H = {et)~^. Moreover, we have from the Einstein Equations that </> ^ \hr^ (t~t ) '^^'^ 
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^ In {t — t^), which resuhs in F = 1 ^(In {t — t*))^, where (J)q is some constant. From this we can have: 

$ c I —, 1^_ y,^o\n{t^U) ^ 



e{t - UY+i[l - ^(ln(i - i,))2] 3a,{t - Uy+i[l - ^(ln(i - U))^]^ 



2/1 8(M^-{l + e)\n{t-U)] \ 

«2 V 1 + e 3(1 + e)3[i _ 4^(in (t - U))^] ) 
8eC0oln(t-t,) / 1 8ee0o[£ +(! + £) In (t-t,)] ^^ 

3[l-^(lnt)2]Vl + e 3(l + e)3[l-4^(ln(i-t,))2]^^ " ' ^ ^ 



St_ , 26/ 1 8e^o[e-(l + 6)hi(t-t,)] 

aie(t-t,)i+^[l-^(ln(i-i,))2] n^\l + e 3(1 + e)3[i _ M^Q^ (^ _ i^))2] 






(48) 



ai(i-i,)i+*[l-^(ln(t-t,))2] 

for long-wavelength case, where the superscript " + " stands for the deflationary region which conies later. 

3. bouncing phase 

When the universe is dominated by the auxiliary scalar 0, the equation of state drops down to below —1, and 
the total energy density of the unverse turns to decrease. When it goes to zero a bounce will happen. It is rather 
complicated to solve the equation (I26p directly, however, we can make some modeling of the evolution to have it 



simplified without losing fidelity. Generalizing the parametrization method in ll| (see also [28|), we parameterize the 
Hubble parameter near the bouncing point of the form; 

H^ait + (3it^ + ... , (49) 

with positive constants ai and /3i of proper dimensions whose magnitudes are determined by the microphysics of the 
bounce. Note that it is unnecessary to contain terms with even power-laws of t. In our case, it can be estimated that 
ai '^ m? . From the relations of variables between Jordan and Einstein frames, we can also obtain the approximate 
value of H with respect to t up to first order of t: 

H = at + higher order terms , (50) 

where a ~ 4ai/9. Since this parametrization is only valid during bounce phase where we can neglect the higher order 
terms of \ri — tjeI with tjb being the comoving time at the bounce point, we obtain the equation for $ of this phase as: 

l>" + 2y(77~?7B)$' + (fc2 + ^)$~0, (51) 

where y = 12aa^/Tr and qb is the scale factor at the bounce point. The solution is: 

where Hi and iFi denotes the Z-th Hermite polynomial and confluent hypergeometric function respectively, with 

3 2y 
The short- and long-wavelength limits of the solution are quite different. For short wavelength (fc2 ^ y)^ it reads: 

, _ .2 fp 
$ ^ (,-" " 2" {^=sin[fc(?7- ?7b)] +Fcos[fc(?7-7/B)]} , (54) 

v2Z 



while for long wavelength (fc2 <C y), it is: 

4 = F + E^{t^ - i^b) - (1 + l)F{7^ ~ ijB? + 0{{v - VB^) , (55) 

where E = -2^+'\/^E/r{-^) and F = F + 2'^E/T{^). 
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4. expanding phase 

The expanding phase can be viewed as time reversal process of the contracting phase. In this period, both the two 
fields will roll down along their potentials and begin to oscillate with a decaying amplitude and redshifted energy 
density. Similar to the case of contracting phase, at the beginning of expansion of the universe, the non-minimal 
coupling term still remains large and dominate over the mass term of the (j) field, and thus (j) can be approximated as a 
slow- rolling field. This drives the total equation of state w to approach —1, and lead to a short period of infiation. But 
soon, when the non-minimal coupling term becomes less important, both fields will oscillate around their minimum, 
behaving like non-relativistic matter again. 

One can also get the short and long wavelength solutions of the metric and field perturbations $, ^ and 6(j) for the 
expanding phase. We only care about the long-wavelength case which make sense for observation today. As similar 
to the contracting case just except for replacing the scripts, one gets: 

^ ^ / 1 8^00 In (^-tt) \^_ 



+ 



+ 



2/ r 1 



■o[e-(l-He)ln(t-tt)] 



+ e 3(l + e)3[l-4^(ln(i-tt))2] 
8ee0oln(i-it) / 1 8£e0o[e+(l + e)ln(t-tt)] X^^- .^q. 



3[1 _ 4|^(ln (i - i^))2] V 1 + e 3(1 + e)3[i _ 4^(in (t _ t^)y 



^^ Si 26/ 1 8g0o[e-(l + e)ln(t-tt)] 



a2(i-it)^+^[l-^(ln(i-it)) 

2a2e{t-^y+i f.^ , ,2ri 4^00,, ,, ,,,21 S^e^Z-o 



1*^^ K^l + e)^ 

for the short infiationary period, and 



((H-.f[l-^(ln(i-tt))^]-^[. + (l + .)ln(t-tt)])^e- (58) 



5^2 ' 




. 61?+ 
' + 5.2 \ 


1 3k^ 



25+ 6D^ 

for the matter-like expanding phase. Here the subscript "e" denotes expanding phase and superscripts " — " and " + " 
stand for inflationary and matter-like regions separately. 

5. spectrum 

Having in hand the solutions in the above sections which stand for different phases, now it is time to connect all of 
them using matching conditions. According to [7] and [8j, we can require that for each point which joins two phases 
together, the three-metric as well as its extrinsic curvature should be continuous. In conformal Newtonian gauge 
which is used in this paper, this indicates that 



* 



V2*-3^(*'-HH$) 

^ 1 



3-HCH' - n^ 



, (61) 



where ± means the difference before and after the transition point. Substituting solutions for perturbations of each 
period into (|6T|) we can get the final results for them. From the last paragraph we can see that at last the solution is 
divided by two parts, one is constant (13+ mode) and the other is decaying (5+ mode). We are only interested in the 
first mode, which dominates over the other one. Since the calculation is rather straightforward and tedious, we only 
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list the final result as follows: 

+ :^[(^-((i)M+(o)Ar+ _ (0)M+(i)iV+ - (i)G+(o)a/+ - (o)G+(i)Af+) 
+ i?-((i)iV+(")iV+ - (i)if+(0)A//+ - (°)if+(i)A/+ + (0)7V+(i)7V+)] 

+ O(^) ■ (62) 

For all the coefficients that appear above, we refer the readers to Appendix B. These coefficients are all the specific 
value at the joint point, so they are independent on k. Moreover, the mode in contracting phase D~ and S~ can be 
easily read from the initial conditions. If we set up with the Bunch-Davies vacuum which implies ^ ~ k~'st~^, we can 
obtain that D~ ^ k^ and S~ ^ k~^. Substituting them into (HU, we can see that at large scales D~ vanishes. The 
contributions of 5*^ in the first order in k^ will be blue-tilted by exactly the right amount to yield a scale invariant 
spectrum. This is the same behavior as the normal case of matter bounce without non-minimal coupling. 

However, different from the normal case, we find that the S~ mode also has contributions in the zeroth order in k"^ 
to the final spectrum. As we know, in the standard case of minimal coupling, the constant mode after the transfer 
point cannot be inherited from the running mode before the transfer. Actually, one may find that the coefficients 
^^'M^ and ^^' M^ vanish except for those terms containing ^ (information from non-minimal coupling term), which 
is consistent with its minimal coupling limit ^ — ^ 0, but when non- minimal coupling terms are introduced, a mixing 
between the two will happen. In our case, since the remaining terms cannot be very large due to the small ^, the 
amplitude of S^ will be suppressed, and in a considerable region of large scales, we can also obtain a scale-invariant 
spectrum [36|. Nevertheless, for extremely large scales, the contribution of S~ of zeroth order in k will become 
non-negligible, and the spectrum will have a red-tilt. 

6. numerical results 

In order to support our analysis above, we also perform the numerical calculation for the perturbations. Fig. [3] is 
the numerical results for the dependence of the metric power spectrum T'^ on the cosmic time t for different comoving 
k modes, where we use the normal definition of power spectrum as: 

'P.{k)^^\M'^^\Dtr (63) 

We choose the zero point on the horizontal axis to be the bouncing point, and set the initial conditions to be Bunch- 
Davies vacuum. One can see from our plot that before the bounce point, the perturbation is dominated by their 
growing modes, and after the bounce, it is dominated by the constant modes, which fits the analytical results very 
well. As for the k dependence, we see that for medium k modes, the power spectrum takes on scale invariance while 
in the extremely small k modes where k < 10~^mpi, the spectrum will present a slightly red tilt. This is because the 
non trivial inheritance of the growing mode in contracting phase to the constant mode in expanding phase at zeroth 
order of k due to the non-minimal coupling effects. Since the scale variance happens only in extreme large scales, we 
expect that it could be tested in the future observational data. 

Furthermore, we calculate numerically the evolution of the particle production S(j) with respect to t and plot the 
result in Fig. |4l Particle production is expected in the region where the squared value of momentum scales of the 
perturbation modes are larger than the potential so that the WKB approximation becomes valid, such as the reheating 
process at the end of inflation z9\. Here at the bouncing point crossing, there is also possibilities to produce particles. 
In the usual case, we can define another variable (p = a^Scj), which satisfies equation of motion as: 

^ + ^lip=..., (64) 

where uf = ^TTm + "^^ ~ i^~ \^^ ■ From above we can see that, although the first two terms are positive-determined, 
there will remain some additional terms due to the geometry of the universe, which would cause w^ to be less than 
zero. If this is the case, the field will have tachyonic behavior and will blow up with particle production (which can 
be called as tachyonic reheating [30|). However, if there is a non- minimal coupling term with a negative coupling 

of motion which becomes: 



constant ^, this effect will get larger since uff, gets a more negative value. This can be seen from the modified equation 



'^+[^+"''~^^~i^'"^^]^^- (^^^ 
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FIG. 3: (Color online) Results of the numerical evolutions of the curvature perturbations with different comoving wavenumbers 
k in non-minimal coupling matter bounce as function of cosmic time t. The initial values of the background parameters are 
the same as in Figure 1. The units of the time axis are M^J,, the comoving wavenumber k is unity for k — nipi. 

This is usually called geometric reheating [3i|. In Fig. |4]we can see that indeed the particle production gets much 
enhancement in the non-minimal coupling case than its minimal coupling counterparts. This is another interesting 
result in our case which are expected to be tested by the future experiments. 
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FIG. 4: (Color online) Results of the numerical evolutions of the field particle production with and without non-minimal 
coupling term for matter bounce. The initial values of the background parameters are the same as in Figure 1. 



IV. CONCLUSIONS AND DISCUSSIONS 



Non-minimal coupling is a very popular subject in cosmology and has been widely studied in the literature. In this 
paper, we have investigated the possibility of generating a matter bounce by the Lee- Wick type scalar field with a 
non-minimal coupling term involved, and studied its perturbations. As in the previous work done in [ll|, a bounce was 
obtained when the non-minimal coupling is not too large to have unexpected effects. However, as it gives correction 
to the effective mass of the field, a short period of deflation/inflation before/after the bounce will happen. 

Using the standard techniques of calculating perturbations of bounce developed from previous works, we have 
calculated the perturbations of the model in detail. One of the big differences from previous one for minimal coupling 
is that the non-minimal coupling to gravity causes the difference between the two scalar perturbations $ and ^, which 
will cause very interesting consequences. We have obtained the solution for each stage of evolution and compiled them 
together using proper matching conditions. We have found that the final dominant mode can be inherited nontrivially 
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from the subdominant mode in the zero-th order of wave number k due to this difference. This will lead to the red-tilt 
of the power spectrum at very large scales, which can be justified in future observation ^^. 

Another result we obtained is that due to the non-minimal coupling, the particle production of the scalar will get 
greatly enhanced. In the region of validity of the WKB approximation when the quantum effects become significant, 
the field will get excited and produce particles. This is very common in reheating process. When there is a negative 
effective mass squared term in the field equation, the particle production will get enhanced due to tachyonic resonance. 
In usual cases, the averaged value of perturbations of the scalar field will get an enhancement through bounce time 
due to geometrical effects. When a non-minimal coupling term is involved in the field equation, the mass of the field 
will be corrected. In the case of negative coupling coefficient, the field will get a more negative effective mass squared, 
and the particle production will get more efficiently enlarged. This is also very interesting phenomenon and can be 
tested by experiments on particle physics. 

Finally we would like to close with some remarks on future works. The bounce with non-minimal coupling is very 
interesting topic since in the early epoch, Einstein's gravity is likely to be modified, which can give some valuable 
effects to the evolution of our universe. It is deserved to pay attention to this project both from the theoretical 
side such as the perturbation, non-Gaussianity, effects on CMB and so on, and from observational side, such as the 
constraints, confirmation or even exclusion from experimental data. All these works are expected to take on in the 
future. 
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Appendix A: Hubble parameter in Jordan and Einstein frames 

This appendix is set to claim the relation between Hubble parameter in the two frames used in the text. As is 
well known, the differentiations of comoving time dr] are the same in the two frames, while that of cosmic time dt are 
not. By definition, we have dt — yFdt where t denotes the cosmic time in Einstein frame. Due to this difference, 
one should be careful when using variables with derivatives in terms of t such as Hubble parameter H. As is already 
mentioned in the previous sections, the Hubble parameter defined in Jordan frame is: 

H^^^^. (Al) 

a adt 

However, we can also define a Hubble parameter in Einstein frame, which is: 

i^-^ = 5, (A2) 

a adt 

where we have introduced the symbol " °" to represent the cosmic time derivative in Einstein frame: °— -^, °°— -^, 
etc. We can see from the definitions that the two Hubble parameters have the same type of relation to the comoving 
Hubble parameters in their frames: H = aH and T-L = ciH . Reminding the relation d — ay/F, it is easy to find the 
description of H in terms of H and F: 

^ ^ djaVF) 
aFdt 

= -^-(^y. (A3) 

y/F WF 

Furthermore, we can derive its derivatives with respect to cosmic time t which are: 

H = —+H( — )--^(^)\ (A4) 

F ^2F' y/F WF 
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Appendix B: coefficients for the perturbation modes 



This appendix indicates all the coefficients that appear in (|62)) . These are calculated using the Hwang- Vishniac 
(Dueruelle-Mukhanov) matching conditions 0, Q : 



(")g+ = 



3aFt3 A 



+ 



B+F[{ 



9aFt 3 aFt 3 A 



-)i? + (16 



4c; 



^ aF^ta aFt3 3aFtt aF^tTi At 



)i/2 



{-SBtAtFH + SA^B^Fff + 6BtCtH^ - 6AtEtH^) 



(Bl) 



(0) t-H- _ 5k^A2 



BtF[^MLH+^2A^+GAV^^ + il. 



5 



n^Aj^l^ 



{-ZB+AtFH + 3AtB+FH + QB+CtH"^ - 6A+E+H^) 



(B2) 



(i)G+ = 



Sc+^( 



5 5-J 

9aFt3 3aFt3 



(-3B^A+FiJ + 3AtB+FH + QB+CtH^ - QA+E+H"^) 



(B3) 



(l)i^^ 



BtF{^ - ^) 



{-3B+A+FH + 3A+B+FH + 6B+C^H^ - 6A+E+H^) 



(B4) 



(o)A/^ 



^+F[(^2v + ^4^)iJ + (16,/i^^ + 



'9aFt^ aFt^A+ 



4e 



4C+ 
■^ aF^fS aFfS 3aFfS aFfS 



)ij2 



{-'SB+AJFH + 'SA+B+FH + 6B+CIH^ - GAJEIH'^) 



(B5) 



(o)7V^ 



.+ i^riMLH. _ (-9 4-^ + fi 4\/3^ - J^ + 36CJ XQ- 



5k^AJF' 



{-3B+A+FH + 3A+B+FH + QB+CtH^ - 6A+E+H^) 



tc 



(B6) 



+ F^ 4e 



(i)M+ 



^^^( 



)aFt? 3aFt3 



{-ZB+A+FH + 3A+B+FH + GB+C^H^ - GA+E+H^) 



(B7) 



(i)7V+ = 






{-3BJA+FH + iAZBJFH + QBtCJH'^ - QAJEJH'^) 



(B8) 



(0) 



M+ = 



K^{2AeFt3 - AC-m^ + ft^eA-F + 2A-HFn - ^eA-RFn + %A- n^S^cfyHneJ ^) 



eim^+UK-^^cPHtiJ^) 



(B9) 



(°)Ar+ = 



K^{2BeFt3 - AE-Rn + ft^eB-F + 2B-HFU - jeB-RFt^ + 8B- K^^(j)RtieJ ^^ 



ei^Ftl + 16K^^(l>Htl 



(BIO) 



(i)Af+ 



«2Fif(-| + |6)A, 



Rei^Fti + IGK^icjyRtiJ^) 



(Bll) 



«7V+ = 



^2^tl(_|+4,)B- 



ReilFti + 16K2^0m§,/3^) 



(B12) 
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where 

" aie(i-t,)i+T[l-M^(ln(t-t,))2] ' " ^2 V 1 + e 3(1 + e)3[i _ M^(ln (i _ t,))2] ^ ' ^ ^ 

C+ = 1 .c^ , (B14) 

ai(t-i,)i+T[l-4^(ln(i-i,))2] ' 

^^_ 86^^{(1 + e)2[l - iifr (In (t-^,))2]-M^[, + (! + ,) In (t-^,)]} 
73^(1 + e)3[l - ^(In (t - t,))2] 



A: 



2e ( 1 8e0o[e-(l + e)ln(t-tt)] 



, 5.7 = 4 (^ + »^'^ol^-U+;i^M^-^t)] A ^ (B16) 



a2t{t-^Y+->[l-^{\u{t-t^)Y\' " n^\l + e 3(1 + e)3[l - ^(In (t - tt)) 



4 ,^2 



n^i(j) 



^' a2{t-t^f+^l-'-^{\n{t-t^W] ' ^^^^^ 



8ef</){(l + e)2[l-lf^(ln(i-it))^]-M^[e+(l + e)ln(i-it)]} 



/3K2(l + e)3[i_i^(ln(t-tt))2] 



(BIS) 



The constants i* and tf are set to guarantee the continuity of H at the joint point tc and te, and one can straightfor- 
wardly check that in the hmit of ^ = 0, we have *^°)M+ and ^°''Af+=0, and a scale- invariant power spectrum can be 
obtained for the whole region of large scales. 
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